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Slender body theory has enjoyed great success in aerodynamics , and 
in recent years there has been considerable interest in applications in 
naval hydrodynamics. Ship hulls are long and slender, and this is a 
logical feature to exploit in simplification of the governing equations. 
Michelles classical theory of wave-making resistance published in 1898 
considered the ship to be a **thin" body and represented the hull by use a 
of distribution of "Havelock sources" over a vertical mid-plane. How- 
ever, ships tend to be "slender" as opposed to "thin" because the beam 
and draft tend to be of the same order of magnitude and both tend to 
be small in comparison to the length. It was hoped, therefore, that 
slender body theory might represent a significantly better alternative 
to Michell’s theory in the area of wave resistance and at the same time 
yield useful results for the hydrodynamic forces in the transverse 
direction both in steady and unsteady motion. Newman (1970) has given 
a review of the application of slender body theory to ship hulls which 
discusses various kinds of hull motion and wave resistance. 

Major advances in the application of slender body theory in ship 
hydrodynamics have been made largely in the past decade or so 
principally by Tuck (1963,1964), Newman (1964,1965), Ogilvie (1969), 
Ursell(1968) and Faltinsen(1971) . The method of approach involves 
the use of matched asymptotic expansions which is a powerful method 
for dealing with singular perturbation problems. Potential flow 
past a slender body is recognized as a singular perturbation problem 
and in such cases an inner and outer expansion is developed. The 
inner expansion is valid in the near-field, i.e., near the body and 
the outer expansion is valid in the far-field. The unknowns re- 
maining in either of the expansion are determined by a matching 
procedure. The general method of matched asymptotic expansions 
is discussed by Van Dyke(1964). 
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In the present report the application of the slender body theory 
to a ship hull in steady motion at high-Froude number is considered. 

While most displacement ships operate in the zero to intermediate Froude 
number range, there is a new generation of high speed ships generally of 
the catamaran configuration which utilize rather thin or slender hulls 
and operate at rather high Froude number. The surface effects ship hulls 
represent the principal example of this type of hull even though the flow 
about an SES hull is somewhat more complex than the one considered here 
because of the pressure difference between the two sides of the hull. 

In the following a single slender hull operating a high-Froude number 
is considered from the view point of slender body theory. 



2. Boundary-value Problem For Ship Hydrodynamics 



We assume the fluid to be ideal (inviscid, incompressible and 
homogeneous) and the motion to be irrotational. If (x,y,z) represents 
Cartesian coordinates fixed in space with the gravitational force 
acting in the negative y-direction and t denotes time, a velocity 
potential (j)(x,y,z,t) can be introduced such that the fluid velocity 
vector is given by 



The condition of incompressible flow applied to (1), i.e.,V.q = 0 
yields Laplace’s equation as the governing differential equation. 



^xx yy zz 
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The fluid pressure P(x,y,z,t) is expressed in terms of ^ and the 
elevation y^ as 



p = ^pV(f). V4)-pgy + C(t) 



(3) 



( 2 ) 



where C(t) may be, in general, a function of time, p denotes the mass 
density of the fluid, and g, the acceleration of gravity! The function 
C(t) is at most a function of time, and can be determined by evaluation 

/s As 

of P+p(J)^+%)Vc()*V({)“fpgy at some reference point. This point is generally 
taken at the free surface at infinity where all of the values are known. 



The surface of the hull is described by the equation s (x,y,z, t)=0. 
Since this surface is considered to be impermeable the kinematic boun- 
dary condition, which imposes the condition that the relative fluid velo- 
city in a direction normal to the hull must vanish, is given by 



^ + V(|)-7s=0 on s=0 (4) 

Dt o t 

On the free surface there are two boundary conditions which must 
be applied. If the elevation of the free surface is described by 
y=n(x,z,t) the kinematic boundary condition obtained in a fashion 
similar to (4) results from setting the substantial deriative of the 
surface n-y=0 to zero, 

]) ^ ^ ^ 

— = 0 on y=n (5) 

In addition to (5) a dynamic boundary condition must also be imposed. 
This boundary condition results from simply applying P=0 on the free 
surface and is given by 

(f)^ + %7({)*7(j)4*gy=C(t) on y=n (6) 

where C(t) is evaluated with the aid of a known reference point in the 
fluid. 



( 3 ) 



In addition to the boundary conditions (4-6) the potential function 
must satisfy certain radiation conditions at x^+y^ <» appropriate to 

the particular problem under consideration. Also, on the bottom the 
kinematic boundary condition 



d) = 0 on y=-h 

y 



is applied if the water depth is finite and denoted by h. 



3. Perturbation Procedure For a Slender Body 



(7) 



At this juncture let us consider a particular problem: a slender 
hull is positioned with its long axis coincident with the x-axis and 
held fixed in space as indicated in Figure 1. A uniform stream of 




Figure 1 Definitions 

magnitude U makes an angle a with respect, to the x-axis so that the 
resulting x,y, and z components are -U cos a, 0, -U sin a, respectively. 
The coordinate system is placed at the mean free surface so that the 
y=0 plane represents the mean free surf ace, and the water depth is taken 
as infinite. The motion will be restricted to steady motion, and the 
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hull is described by s(x,y,z;e) = 0 where the parameter e = beam or 
draft -to -length ratio. 

Outer Problem 



The problem which is here defined is recognized to be dependent 
on two dimensionless parameters, the beam/length ratio, e, and the 
incidence angle, a. Thus, as an initial step we may expect the solu- 
tion to be expandable in a power series in both e and a. (At least 
the solution should behave in a regular manner with the exception of 
the region near the body.) Introducing the body potential, the 
total potential may be expressed as 

(j) = (j)Cx,y,z;e,a)-Ux cos a-Uz sin a (8) 

where the body potential <J) is expressed by 

(J>=e(f>^+a(j)%0(£2)+0{a2) (9) 

in the outer region. The first term of the Taylor series, (0,0)=<i>Q, is, 
of course, zero because if both e=0 and a=0 the uniform stream is 
directed parallel to a line of zero thickness and the body potential 
vanishes. Thus, expanding sina and cosa about a=0 in (8) and substi- 
tuting into (9) gives the total potential as 

()) = £())^ + a((J)^-Uz)-Ux+0(£^)-K> (a^) (10) 

and substituting this form of the potential into (2) shows that 

V2(j)^=v2^“ =0 (11) 

since a and e are independent parameters. 
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In similar fashion, we may express the free surface elevation in 
the form of (9) , 



n=eri^-^ari^+0( e^)+0(a^)40(ae) (12) 

Substituting (8), (9) and (12) into (5) and collecting coefficients of 
like powers of e and a gives: 



and 

^ 

The dynamic free surface boundary condition is obtained by use of 
(6) by evaluating C(t)=h\]^ at x=«> on y=0 and substituting in the form 
of <() and n as given in (12) and (8). Here again two separate boundary 
conditions result because of the independence of the two small parameters, 
e and a. The first of these is obtained by setting a=0 as 






iL 

g 



= 0 on y=0 



(15) 



and the second by setting e=0, 



JL.^" = 0 

1 g ^Ix 



Eqs. (13) and (14) together with 
surface boundary conditions, 



^1 

ly g ixx 



on y=0 (16) 

(15) and (16) yield the linear free 



= 0 



(17) 



and 




U2 . a 
g 



0 



(18) 
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Eqs. (17) and (18) represent familiar linear boundary conditions for 
the free surface. If we now let the Froude number U^/gL (where L is 
considered to be 0(1)) become large, (17) and (18) become simply 



((>1 = 0 



on y=0 



(19) 



(j)^ = 0 on y=0 (20) 

In summary, we have derived a boundary-value problem for the outer 

potential. The solution is singular in the region of the body so the 

hull kinematic boundary condition is not applied to the outer solution. 

0 

We, therefore, have in summary for (J)^ , 



£ £ £ ^ 
(j). + <(>T ((>T ” 0 

ixx lyy izz 



<}), = 0 on y=0 



and similarly for (j)^. 



((>? + + (j>? = 0 

-XX lyy Izz 



d), = 0 on y=0 



( 22 ) 

(23) 



(24) 



(25) 



The boundary condition on the hull is in this case to be replaced by a 
matching condition. 



Inner Problem 



We now consider the inner potential which is valid near the slender 
body. The solution to the slender body problem is singular in the inner 
region but the degree of singularness at the outset is unknown. We there- 
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fore begin by writing the body potential as 



(}> = Aq (e)$Q (x,Y,Z;a)+Aj (£)$J (x,Y,Z, ;a)+. . . 



( 26 ) 



where the Y and Z coordinates are stretched by the parameter e 



Y=y/e 



(27) 



Z=z/e 



(28) 



and (e) , (e) , etc. are called gage functions as defined by Van Dyke 
(1964). These gage functions will eventually be determined by use of 
the principle of least degeneracy also discussed by Van Dyke. 

The body potential is expressed in (26) as an inner perturbation 
series in e, the beam-length ratio. Here the incidence angle, a, is 
allowed to remain finite and is not considered in any limiting process. 
We shall return to this point later. 

If the inner representation of the body potential given by (26) 
is substituted into (8) and the result substituted into the Laplacian, 
(2), we find, 



Clearly, if Ai(e)/A 2 (€) is at least of 0(c) then the three dimen- 
sional Laplacian degenerates to the two-dimensional form for the leading 
term in the inner expansion: 



The inner representation, being valid near the hull , must satisfy the 
boundary condition on the hull. If s is expressed in inner variables as 




$ + $ =0 



(30) 




( 8 ) 



s(x,y,z;e) = S(x,Y,Z;e) = 0 then the hull kinematic boundary condition, 
which is obtained by use of (26) and (8) with (4), gives 



(■ 



X 






)• 



(g) 

A^(e) 






S +$i S +$ 



X 




E^U eUaS„ 

X — ij 

A (e) 

o o 



(31) 



Now if a were to be considered to be fixed at 0(1) then the principle of 

least degeneracy would require that A (e) = e, and the hull kinematic 

o 

condition would read 



^ vS =UaS„ on S=0 

oY Y oZ Z Z 



(32) 



This boundary condition simply indicates that in any transverse plane 
perpendicular to the x-axis the body potential must just cancel the 
cross-flow velocity, Ua. It is noted that in such a case the cross- 
flow is an order of magnitude larger than the perturbation associated 
with the thickness distribution of the body and thickness effect is, 
therefore, lost. 

Alternately, if the angle of attack were taken to be zero then the 
principle of least degeneracy leads one to A^(e) = e^ and (31) gives, 

^ S S =US on S=0 (33) 

OY Y OZ Z X 

a familiar result in slender body theory. 

A third possibility is to take a to be of the same order of magni- 
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tude as e and retain both of the terms on the right hand side of (32) 
and (33). This gives A^(£)=e^ and 



An alternate way of writing (34) is 



= n/'Sx . 

^ VSn e oZj 



(35) 



where N=^/e and n represents distance measured normal to the curve in the 
Y^Z plane described by the equation S(x,Y,Z;e) = 0 with held fixed. 



The free surface condition to be applied to the inner solution is 
obtained by eliminating q between (5) and (6) • The boundary condition 
obtained in this way is 

(> (V<i>*V(|) )+(f> (V(Ji*V(fi )+g(J) =0 on Y=n/£ (36) 

X X 2> z Y 

Substitution of (26) and (8) into (36) gives 






(37) 



Now 



, if e-H) with a=fixed, (37) gives 



$ _ =0 on Y=0 
oZZ 



(38) 



or from this we can deduce that 



$ =0 on Y=0 
o 



(39) 



is the appropriate boundary condition in the inner region. The problem 
for the first term of the inner expansion for the high-Froude number 
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slender body theory is, in summary, 



oYY oZZ 



(40a) 




TT/ X 3N . ^ ^ 



(40b) 



<I>Q =0 on Y=0 (40c) 

There is no radiation condition specified in (40) ; this condition is 
replaced by matching with the outer solution. 



4. Solution for the Inner Problem 

The inner problem specified in (40) is essentially a two-dimensional 
problem as indicated in Figure 2. Such problems can be solved for rather 




Figure 2 Inner Problem 

arbitrary shapes by use of a Greenes function of the form 



G(Z,Y;C,n)=ln>/(Z-?)2+(Y-n)2-ln\/(Z-02 +(Y+n)2 (41) 
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The first term in (41) represents a source, and the second term represents 
a negative image with respect to the Z axis. The potential 4>^may be 
represented by a surface distribution over the hull of these two-dimen- 
sional singularities. 



In particular, the potential may be represented by 



$ = 
0 



2ir y^F(x;ri,i 



,?)G(Z,Y;?,n)dC (42) 

where C denotes the curve in the Y-Z plane representing the cross- 
section of the hull and dC denotes an arc length along that curve. 
Using (41) for G, (40a) and (40c) are automatically satisfied. The 
hull boundary condition (40b) then yields the integral equation for 
the source strength distribution, F(x;Y,Z) as 



F(x;Y,Z)+^ /F(x;n,?) (Z,Y; 4 ,n)dC 

II f Q oIN 

\% = 



on C 



(43) 



This equation can be solved numerically by subdiving C into a finite 
number (say N) of segments. Eq. (43) can then be reduced to a set of 



N linear equations in N unknown values of over each of the segments. 



Cylindrical Hull 

For the case of a semi-immersed circular cylinder, a multipole 
series may be used to represent the solution. Consider the surface 
described in cylindrical coordinates by 



S(R,0,x)=R-R^(x) 



(44) 



where R-r/e denotes the radial coordinate in inner variables and 
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denotes the radius of the hull in the inner variable (R^ 
is, in general, a function of x. The kinematic boundary condition on 
the hull, (40b), takes the form, 

d) 1 

= U(R^ + I cose) on R=Rj^ (45) 

where the polar coordinate system is related to the cartesian system by 



Z=R cos 0 (46a) 

Y=-R sin 0 (46b) 



The solution for $ satisfying both (40a) and (40b) is given by 

0 

CO 

$ =5^ sin(n0) 

n=l (47) 

Application of the hull boundary condition (45) now gives 
00 a n ^ 

E :(n+l) ne=-UR^ - T U cos0 (48) 

n-l 



(49) 

if m?^n (50) 

(51) 

(52) 

we find 



By use of the results 



J sin(n0)d0-^ cos^(n 0 )d 0 = 7 r /2 
0 o 

^TT 

y sin (n0 ) sin (m0 ) =/ cos (n0 ) cos (m0 ) d0=O 



/ 

/ 

o 

7T 



/ 



2n 

sin(n0)cos0 d0= — 7 —, 5 
n"^-l 



si„(„9)de- 



n 



n 



n 
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__a 
£ n 






COS0 sin(n0)d0, n=2,4,6,,... (^3) 



or, carrying out the integration, 






and 



"-T (:?^) ■ 






(54) 



= _ 2 



n 



1 > 3 , 5 , • • « ( 55 ) 



The inner potential function is, therefore, represented by the series 



^ o_ ^ 
0 e 



(ii) 



. . .. 

^n— sin(ne) 



n=2,4 



(56) 



n=l,3,.. ' 



T. ^+1 . / AX 

sin(n0) 



The dynamic pressure on the hull can be obtained from Bernoulli's 
equation, disregarding terms of O(e^) and O(e^ct^), 

P-K^Joto^- (»2^ i^S-(|)M*^^si„e+(a)%2sin29j (57) 



The normal (or yaw) force per unit length along the hull can 
be determined through integration of the pressure over the immersed 
surface. That is, 

7T 

(58) 



dF^= P Ry COS0 d0j dx 



( 14 ) 



After substitution of (56) into (58 ) and some lengthy manipulations 
and integrations we get, for the normal force on the segment dx. 



dF^= -£2 qpU2 ^ Rh Rjj* 2+^ (izlzlf: n(l+(-l)^ 

n=2 I ^ 



>dx (59) 



The net normal force is obtained by integration from the bow where 
Rjj=0 to the stern where 



-a ^ pu^ ^ -r: 

^ L 2 “S ^ 



“i: 



n(l+(-l)^) 
n^-l 



(60) 



n=2 



The dimensionless normal (or yaw) force coefficient based on the stern 
area, therefore, becomes 
-F^ 8 



f = 



2 Jjpu^TT r2 ^ 
— lie 
2 S 



2+ 



E 

n=2 



l-(-l)'' n(l+(-l)^) 



n 



u2-l 



(61) 



where r„ is defined as R„ =r„ / £ and ru represents the hull radius in outer 

Hg -Hg Hg 

variables. For boat-tailed hulls we get the familiar result F2=0 because 

Rhs=o- 

The summation in (61) can be somewhat rearranged to read 

CO 

(62) 



TT^I 3 ^ j n(n+2) 

n=3,5,7 



The series converges rapidly to approximately 0*128 so 



f =0.75a 
z 



(63) 



The value of f^=0.75a for the cylindrical hull at high Froude number can 
be compared with the classical value of f^=2a for the same hull operating 
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with the low-speed free surface boundary condition, D4>^/3Y=0 on Y=0. 

It is reasonable to expect a somewhat smaller value for f^ for the 
high-Froude number case because the high“Froude number free surface 
boundary condition makes the Y=0 plane act like a zero pressure sur- 
face and, therefore, right at the free surface there can be no pres- 
sure difference across the hull. The low Froude number (actually low 
to moderate) boundary condition makes the free surface act like a rigid 
wall and a pressure difference across the hull can exist at the Y=0 
level. 



The moment about the Y-axis is given by 

b 



dM = -/ xdF 

y J 2 



x=-a 



(64) 



where a and b denote the distance from the origin to the stern and bow, 
respectively, so that the hull length is L=a+b. Thus, using (59) with 
(64) gives 



^ = 



2 





n(l+(-l) 






(65) 



where 



Ag=area of the stern section below the mean waterline 
¥=hull displacement 
L=hull length 



When the series is evaluated we get 



my =0 . 75 




( 66 ) 
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and for boat-tailed (or closed) bodies the moment coefficient becomes 
my=0,75a. This value may be contrasted with the low Froude number 
result of 2.0a which is often referred to as Munk’s moment. 

5. Conclusions 

Slender body theory has been applied to a slender ship hull opera- 
ting at high-Froude number. For the case of a semi-circular hull, closed- 
form results were obtained for the yaw force and moment. These values 
are less than half the values of the equivalent coefficients associated 
with the rigid-wall free surface boundary condition. 
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